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Abstract 

Thin  (~30-45  jim)  gadolinia-doped  ceria  (GDC)  electrolyte,  Ni  +  yttria- stabilized  zirconia  (YSZ)  anode- supported  cells  were  fabricated  with 
platinum  reference  electrodes  (probes)  embedded  within  the  electrolyte  layer.  The  exposed  parts  of  the  probes  were  sealed  off  from  the  atmosphere 
using  an  oxygen  impermeable  glass.  Thus,  no  oxygen  exchange  could  occur  and  voltage  measurements  at  the  probes  correspond  to  the  reduced 
negative  electrochemical  potential  of  electrons  (<p  =  —pie/e,  where  /xe  is  the  electrochemical  potential  of  electrons  and  e  is  the  electronic  charge). 
The  <p  was  measured  at  650  °C  under  open  circuit  conditions  with  air  circulated  past  the  cathode  and  hydrogen  circulated  past  the  anode.  The 
open  circuit  voltage  was  ~0.84  V.  The  voltages  between  the  anode  and  reference  electrodes  were  very  small — less  than  30  mV.  This  was  attributed 
to  significant  electronic  conduction  throughout  the  electrolyte  from  anode/electrolyte  interface  to  cathode/electrolyte  interface.  Thus,  voltage 
drop  mainly  occurs  across  the  cathode/electrolyte  interface.  The  variation  of  chemical  potential  of  oxygen,  /xq2,  within  the  GDC  electrolyte  was 
estimated,  which  showed  a  large  drop  in  /x02  across  the  cathode/electrolyte  interface.  The  oxygen  partial  pressure  in  GDC  just  inside  the  cathode 
was  estimated  to  be  ~10-22  atm.  Thus,  in  GDC-based  anode- supported  cells  the  effective  electrolyte  thickness  is  very  small. 

©  2008  Published  by  Elsevier  B.V. 
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1.  Introduction 

1.1.  Transport ,  local  equilibrium  and  ionic  and  electronic 
currents 

Numerous  studies  on  transport  through  ionic  conductors  have 
been  reported  over  the  years.  Oxygen  ion  conductors  have 
been  of  particular  interest  due  to  their  applications  in  such 
devices  as  solid  oxide  fuel  cells  (SOFCs),  oxygen  separators,  and 
electrolyzers.  The  materials  studied  include  those  exhibiting  pre¬ 
dominantly  ionic  conduction  (transference  number  >0.99)  and 
those  exhibiting  significant  electronic  transport  in  addition  to 
ion  transport  (mixed  ionic  and  electronic  conduction  (MIEC)). 
Yttria-stabilized  zirconia  (YSZ)  is  essentially  a  pure  oxygen  ion 
conductor,  while  rare  earth  oxide-doped  ceria  exhibits  mixed 
ionic  and  electronic  conduction.  The  basis  of  the  transport  the¬ 
ory  was  developed  over  50  years  ago  [1-3] .  Virtually  all  reported 
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studies  on  solid-state  transport  through  ionic  and  mixed  ionic- 
electronic  conductors  are  based  on  two  assumptions  [4-6] : 

( 1 )  Transport  of  a  charged  species  occurs  down  its  electrochem¬ 
ical  potential  gradient,  namely 

ji  =  —CiBiV  fli  (1) 

where  Q  is  the  concentration  of  species  i\ 
Bi  =  Dfk b  T  =  mobility  of  species  i,  where  D[  =  diffusivity 
of  species  i\  jli  is  the  electrochemical  potential  of  species  i. 

The  corresponding  current  density  due  to  species  i  is  given 
by 

li  =  “  — V/Z;  (2) 

Zie 

where  07  is  the  conductivity  due  to  species  i  and  Zi  is  the 
valence  of  species  i. 

(2)  Local  equilibrium  exists  everywhere  in  the  system  which 
allows  one  to  write  all  transport  equations  in  terms  of  gradi¬ 
ents  in  chemical  potentials  of  neutral  species  and  electric 
potentials  (electrochemical  potential  of  electrons,  /xe,  or 
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Nomenclature 

B[  mobility  of  species  i  (velocity/unit  force) 

C\  concentration  of  species  i 

D{  diffusivity  of  species  i  (cm2  s-1  or  m2  s-1) 
e  electronic  charge  (Coulombs) 

E  Nemst  voltage  (V) 

Em  open  circuit  voltage  (V) 

E3-  internal  EMF  across  anode/electrolyte  interface 

in  the  equivalent  circuit  (V) 

Ec  internal  EMF  across  cathode/electrolyte  interface 

in  the  equivalent  circuit  (V) 

Eel  internal  EMF  across  electrolyte  (just  inside  inter¬ 
faces)  in  the  equivalent  circuit  (V) 

V  cell  voltage  (V) 

rf  anode/electrolyte  interfacial  ionic  area-specific 

resistance  (Q  cm2) 

rf  cathode/electrolyte  interfacial  ionic  area- specific 

resistance  (Q  cm2) 

rf  electrolyte  area- specific  ionic  resistance  (£2  cm2) 
rf  anode/electrolyte  interfacial  electronic  area- 

specific  resistance  (Q  cm2) 

rf  cathode/electrolyte  interfacial  electronic  area- 

specific  resistance  (Q  cm2) 

R{  =  rf  +  rf  +  rf  cell  area-specific  ionic  resistance 
(Q  cm2) 

Re  =  rf  +  rf  +  rf  cell  area-specific  electronic  resistance 
(Q  cm2) 

7e  electronic  current  density  through  the  membrane 

(A  cm-2) 

7i  ionic  current  density  through  the  membrane 

(A  cm“2) 

4  measured  current  density  (A  cm-2) 

Boltzmann  constant  (J  deg  - 1 ) 
po2  oxygen  partial  pressure  (atm) 

R  ideal  gas  constant  (J  deg-1  mol-1) 

7?l  load  (£2  cm2) 

Greek  symbols 

li i  chemical  potential  of  species  i  (Jmol-1  or 

J/species) 

jli  electrochemical  potential  of  species  i  (J  mol-1  or 
J/species) 

cp  =  —jfQ/e  reduced  negative  electrochemical  potential  of 
electrons  (V) 


negative  reduced  electrochemical  potential  of  electrons,  (p  = 
—jle/e,  which  is  also  the  reduced  electrochemical  poten¬ 
tial  of  electron  holes).  For  a  cation  with  valence  zi  (zi  >  0), 
namely  MZ1 ,  the  local  equilibrium  is  described  by  [2]: 

MZl  +  z\e  M  (3) 

which  gives 

A mzi  +  z\  Ae  =  Mm  =  Amzi  -  Ziecp  (4) 


where  /xm  is  the  chemical  potential  of  neutral  M.  Similarly, 
for  an  anion  of  valence  zi  (or  —\zi\  since  Z2<0),  namely 
XZ2  or  X-^,  local  equilibrium  is  described  by  [2]: 

X+\z2\e=iX2  +  \z2\eOX-W  (5) 

which  gives 

MX  +  |Z2lAe  =  jMX2  +  IZ2lAe  =  MX-^I 

=  2MX2  +  Z2 e(p  =  2 MX 2  -  \z2\ecp  (6) 

where  fix  is  the  chemical  potential  of  neutral  X  and  fix2  is 
the  chemical  potential  of  neutral  X2. 

These  assumptions  are  universally  made,  either  explicitly  or 
implicitly  in  the  reported  studies  on  solid-state  transport  through 
ionic  and  MIEC  conductors.  Focal  equilibrium  also  implies  that 
any  changes  in  thermodynamic  potentials  from  the  equilibrium 
value  occur  such  that  [7] 

<$Amzi  +  Zl^Ae  —  =  $Amzi  —  Z\e8cp  (7) 

and 

S/Xx  +  \Z2 1  $Ae  —  58/XX2  +  \Z2\  SAe  =  SAx~lz2l 

=  28MX2  -  \z2\e§(p  (8) 

where  8/x  denotes  a  variation  in  the  thermodynamic  potential  fi. 
The  application  of  the  local  equilibrium  concept  to  a  predom¬ 
inantly  ionic  conductor  has  important  implications  concerning 
transport  and  chemical  potentials  (stability)  [7]. 

The  net  current  density  due  to  all  species  is  given  by 

/  =  VAi  (9) 

i  izie 

assuming  the  electrodes  can  exchange  all  species  (not  blocking 
to  any  of  the  species).  If  for  example  the  material  is  a  predom¬ 
inantly  ionic  conductor,  then  current  only  due  to  ionic  species 
dominates.  That  is,  7 = 7ion  +  7eiec  =  Aon  since  |7ion  |  >  |/eiec  I  •  That 
is,  in  such  a  case,  virtually  all  of  the  current  is  due  to  ionic 
species  (typically  one  species)  and  transport  of  the  remaining 
species  can  be  ignored,  insofar  as  the  magnitude  of  the  current 
is  concerned.  This,  for  example  is  the  case  with  YSZ  as  the  elec¬ 
trolyte  in  which  the  dominant  current  is  due  to  oxygen  ions.  It 
is  a  common  practice  to  equate  all  of  the  measured  current  with 
the  ionic  current  through  YSZ  membranes. 

An  important  point  to  note  is  that  electronic  current  in  a 
predominantly  ionic  conductor  can  be  neglected  only  when 
considering  the  total  current.  However,  insofar  as  the  local  equi¬ 
librium  is  concerned,  electronic  current  cannot  be  neglected,  no 
matter  how  small  the  electronic  current  is.  This  is  because,  the 
conditions  described  by  Eqs.  (7)  and  (8)  can  be  satisfied  only 
if  electrons  can  be  transported,  in  and  out,  so  that  variations  in 
all  thermodynamic  potentials  can  occur  under  any  applied  or 
internal  stimuli  [7].  That  is,  the  electronic  conductivity  cannot 
be  set  identically  to  zero  or  the  electronic  resistivity  cannot  be 
set  to  infinity.  As  has  been  shown  previously  this  is  because  the 
electronic  current  enters  the  relations  for  the  chemical  potentials 
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of  neutral  species  as  a  multiplicative  factor  (a  product  of  elec¬ 
tronic  current  and  electronic  resistance)  [7-9].  For  oxygen  ion 
conductors,  this  means  even  in  a  predominantly  ionic  conductor 
such  as  YSZ,  the  electronic  current  cannot  be  neglected  since  it 
enters  as  a  multiplicative  factor  (multiplication  of  the  electronic 
current  and  the  electronic  resistance)  in  the  equations  for  oxy¬ 
gen  chemical  potential,  /xo2  •  Neglect  of  electronic  current  leads 
to  indeterminacy  of  the  chemical  potentials  of  neutral  species 
and  violates  the  fundamental  local  equilibrium  criterion  [7,8]. 
The  preceding  implies  that  while  the  magnitude  of  the  electronic 
current  may  be  small  (as  an  additive  term  in  the  total  current), 
its  effect  on  the  chemical  potential  is  multiplicative  and  can  in 
fact  dictate  the  very  stability  of  the  material  under  the  given 
experimental  or  operating  conditions. 

7.2.  Role  of  interfaces 

The  electrode/electrolyte  interfaces  are  usually  physically 
distinct,  such  as  a  sharp  change  in  micro  structure  across  an 
interface.  Properties  of  the  electrolyte  in  close  proximity  to 
the  interface  are  different  from  those  in  the  bulk  electrolyte 
due  to  compositional  variations  (segregation  of  some  species, 
and  depletion  of  other  species),  space  charge  effects,  etc.  These 
regions  can  extend  several  nanometers  into  the  electrolyte  and 
in  some  cases  even  several  microns.  They  may  or  may  not  be 
present  in  all  materials.  Even  in  the  absence  of  such  gradi¬ 
ents  (compositional,  space  charge),  the  electrolyte  region  very 
close  to  the  interface  must  always  exhibit  different  proper¬ 
ties  from  the  bulk.  This  is  because  atomic  bond  lengths  of 
near  surface  atoms  are  different  from  those  in  the  bulk  due 
to  unsatisfied  bonds  on  the  surface.  For  this  reason,  transport 
properties  of  interfacial  regions  must  also  always  be  differ¬ 
ent  from  those  in  the  bulk.  Since  electronic  conductivity  of 
the  bulk  electrolyte  cannot  be  identically  zero  for  the  afore¬ 
mentioned  local  equilibrium  criterion,  it  means  that  electronic 
conduction  must  also  occur  through  the  interfacial  regions. 
That  is,  interfacial  regions  must  transport  both  ionic  and  elec¬ 
tronic  species,  no  matter  how  small  is  the  current.  Thus,  the 
electrolyte  for  any  modeling  purposes  must  be  divided  into  at 
least  three  separate  regions:  the  two  interfacial  regions  and  the 
bulk. 


When  the  two  surfaces  of  a  membrane  are  exposed  to  two 
different  atmospheres,  such  as  fuel  and  oxidant,  the  chemical 
potentials  vary  through  the  membrane — that  is,  they  are  func¬ 
tions  of  position.  Since  defect  concentrations  are  functions  of 
chemical  potentials,  the  transport  properties  are  also  functions 
of  position  within  the  membrane.  For  this  reason,  the  spatial 
variation  of  chemical  potential  is  expected  to  be  nonlinear.  If 
the  ionic  and  electronic  transport  properties  are  assumed  to  be 
independent  of  position,  the  spatial  variation  of  chemical  poten¬ 
tials  is  expected  to  be  linear.  As  was  discussed  recently  such  an 
assumption  only  changes  the  details — but  not  the  general  fea¬ 
tures  of  transport  [7] .  In  this  manuscript,  this  assumption  is  made 
for  illustration  purposes  only  since  it  allows  the  use  of  simple 
equivalent  circuits,  such  as  the  one  shown  in  Fig.  1.  The  equiv¬ 
alent  circuit  shown  in  Fig.  1  ignores  concentration  polarization 
and  activation  polarization  associated  with  the  extended  reaction 
zones  in  the  porous  electrodes  [7,9].  The  various  parameters  in 
the  equivalent  circuit  are  as  follows:  rf,  rf ,  and  rf  are  ionic  spe¬ 
cific  resistances  of  the  cathode/electrolyte  interface,  electrolyte, 
and  electrolyte/anode  interface,  respectively;  r£,  rf,  and  are 
electronic- specific  resistances  of  the  cathode/electrolyte  inter¬ 
face,  electrolyte,  and  electrolyte/anode  interface,  respectively; 
and 


EC  _  /4  -  /^2 

4e 

Ezl  =  **02  ~  ^02 

\e 

and 

EZ  _  ~  ^02 

\e 


(10) 

(ID 


(12) 


are  the  internal  EMFs  across  the  cathode/electrolyte  interface, 
electrolyte,  and  electrolyte/anode  interface,  respectively,  which 
are  given  in  terms  of  the  various  oxygen  chemical  potentials, 
/Zq2  s.  The  sum  of  the  three  corresponds  to 

Ec  +  £el  +  £a  =  ^2 - =  e  (13) 

4e 


which  is  the  Nernst  voltage. 


Fig.  1.  A  simplified  equivalent  circuit  for  an  SOFC.  The  ri’s  are  the  ionic  resistances,  rf  s  are  the  electronic  resistances,  and  Ea’ s  are  the  internal  EMF’s  (representing 
A/xo2’s).  The  circuit  elements  between  the  two  filled  circles  are  not  physically  separable  and  collectively  represent  the  equivalent  circuit.  Rl  denotes  the  external 
load.  The  open  circuit  condition  is  given  by  /?l  oo. 
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1.3.  Ionic  and  electronic  currents 

For  an  SOFC  with  YSZ  as  the  electrolyte,  the  electronic 
resistance  is  much  greater  than  the  ionic  resistance.  With  ref¬ 
erence  to  the  equivalent  circuit  given  in  Fig.  1,  this  means 
R&  =  r£  +  rf  +  rl  rf  +  rf  +  rf  =  R[.  Thus,  the  magnitude 
of  the  ionic  current  is  much  greater  than  the  magnitude  of  the 
electronic  current,  that  is,  |/i|  |/e|,  and  the  current  measured 

in  the  external  circuit  7l  =  h  +  7e  is  predominantly  due  to  the 
ionic  current  through  the  membrane.  The  general  equations  of 
transport  are  given  by  [1-3]: 

Iq2-  =  h=  V A02-  =  -^-V^02  -  <702-  V<?  (14) 

Iq  =  —<JQ  Wcp  (15) 

In  terms  of  the  various  area  specific  resistances,  the  currents 
are  given  as  follows: 

Cathode/electrolyte  interface  (assuming  a  one-dimensional 
case,  such  that  all  potentials  are  functions  of  distance  perpen¬ 
dicular  to  the  electrolyte  surface): 

/C  = _ j_  /4  ~  lAh 

O2-  4e  rf  rf 

Bulk  electrolyte: 

.el  _ L  ~  11  (h  f 

O2-  “  4e  rfl  rf1 

Anode/electrolyte  interface: 

4  . .  a.  1 1 II  a  II 

/a  =  1  Mq2-Mq2 

O2--  4e  rf  rf  ’ 

In  steady  state,  the  divergence  of  ionic  and  electronic  cur¬ 
rents  is  individually  zero,  that  is  V  •  7Q2-  =  0  and  V  -7e  =  0, 
or  dIQ2-/dx  =  0  and  d7e/dv  =  0  in  a  one-dimensional  case, 
which  means  the  ionic  current  is  constant  (spatially  independent) 
through  all  three  regions,  7^2_  =  lf2_  =  7^2_  =  IQi-  =  7i;the 
electronic  current  is  also  constant  through  the  three  regions, 
If  =  if  =  If  =  7e.  Eqs.  (16)— (18)  can  be  rearranged  to  write 
the  chemical  potentials  of  oxygen  just  across  the  interfaces 


inside  the  electrolyte,  namely  /Zq2  and  /Zq2  in  terms  of  ionic 
and  electronic  currents  [7,9]. 

Mo2  =  Mo2  +  Mrf /,  -  rec/e)  (19) 

l4h  =  Mo2  -  4e(rf/i  -  rea/e)  (20) 

Also,  the  rearrangement  of  equations  for  electronic  transport 
gives 

<PC  =  (P1  -  rceIe  (21) 

(pa  =  <p11  +  rea/e  (22) 


As  a  general  rule,  there  should  be  abrupt  changes  in  both  /zo2 
and  ip  across  both  interfaces  (since  interfacial  regions  are  very 
thin — on  the  order  of  nanometers,  but  with  different  properties 
than  the  bulk).  Fig.  2  shows  schematic  variations  of  /zq2  and  <p 


rc  =V]  -  <PC 


(16) 


7el  =<PC~<pa 


(17) 


.a  _  <P&  ~  A1 

1  P 


(18) 


Interface  Interface 


Fig.  2.  Schematic  variations  of  /xo2  and  (p  across  a  membrane  corresponding  to 
the  equivalent  circuit  in  Fig.  1.  If  the  transport  properties  of  the  membrane  are 
dependent  upon  /xo2  (and  thus  on  position),  the  variations  of  /zo2  and  (p  will  be 
nonlinear  (shown  by  broken  lines). 

through  an  SOFC.  Relative  values  of  transport  properties  (elec¬ 
trolyte  and  interfaces)  can  substantially  influence  the  variations 
of  /xo2  and  p  [7].  Also,  if  the  transport  properties  of  the  elec¬ 
trolyte  are  functions  of  /xo2,  then  the  variations  of  [i o2  and  <p 
through  the  electrolyte  will  be  nonlinear  (shown  by  broken  lines 
in  Fig.  2).  In  Ag2S,  for  example,  the  variations  of  /ZAg  and  cp  are 
nonlinear  [1]. 

1.4.  Review  of  measurement  of  potentials  using  selective 
probes  and  calculations  of  potentials 

Several  studies  have  been  reported  on  the  estimation  of 
thermodynamic  potentials  in  solid-state  ionic  and  mixed  ionic 
electronic  conductors  [8,10-22].  Analysis  of  a  selective  probe 
for  the  measurement  of  cp  and  chemical  potential  of  neutral  sil¬ 
ver  in  Ag2S  was  first  reported  by  Hebb  [1].  Depending  upon 
whether  the  tip  of  the  probe  can  exchange  electrons  or  silver 
ions,  the  measured  voltage  could  be  correlated  with  <p  or  /ZAg. 
After  this  initial  work  by  Hebb,  several  researchers  have  used  the 
selective  probe  method.  For  example,  the  measurements  of  ther¬ 
modynamic  potentials  in  Ag|AgBr|Pt  and  Ag|AgBr|  |AgCl|Ag 
have  been  reported  [19].  For  the  accurate  measurement  of  /XAg 
in  Agl,  Pt  probe  shielded  by  an  alumina  tube  has  been  used 
in  one  study  so  that  only  the  tip  contacted  Agl  [20] .  For  mixed 
ionic  electronic  conducting  CU2O,  MIEC  probes  have  been  used 
[10]. 

In  cation  conductors,  as  the  ambient  atmosphere  does  not 
contain  the  neutral  metal  corresponding  to  the  cation,  metal 
exchange  with  the  atmosphere  is  generally  not  an  issue.  In 
anion  conductors,  especially  oxygen  ion  conductors,  interac¬ 
tion  with  oxygen  in  the  ambient  atmosphere  must  be  taken  into 
account.  The  concept  of  embedded  reference  electrodes  has  been 
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explored  by  Reed  et  al.  [21].  In  their  studies,  a  Pt  probe  was 
embedded  in  an  YSZ  sample  ~200  p,m  in  thickness.  However, 
it  was  not  stated  whether  or  not  the  probe  was  isolated  from 
the  atmosphere.  An  experimental  study  on  the  measurement 
of  oxygen  chemical  potential,  /x o2,  on  a  thick  gadolinia-doped 
ceria  (GDC)  sample  has  been  reported  by  Mineshige  et  al.  [22]. 
A  4-mm  thick  sample  of  GDC  was  exposed  to  fuel  on  one 
side  and  air  on  the  other  side  at  1000  °C.  After  equilibration, 
the  sample  was  sectioned  and  the  /x o2  profile  was  obtained 
using  Raman  spectroscopy  [22].  Data  were  obtained  in  inter¬ 
vals  ranging  between  ~50  p,m  and  ~400  [Jim.  The  interval  was 
~400  Jim  near  the  anode  and  ~50  [xm  near  the  cathode.  Sev¬ 
eral  studies  have  been  reported  on  theoretically  estimating  the 
/xo2  variation  within  solid  electrolytes  [11-16].  However,  in 
all  of  these  studies  it  was  assumed  that  the  /xo2  in  the  solid 
electrolyte  in  contact  with  the  gas  phase  equilibrated  across 
the  interface.  As  has  been  shown  recently,  equilibration  of 
/x o2  across  gas/solid  interface  cannot  be  assumed  when  the 
two  surfaces  of  a  solid  are  exposed  to  atmospheres  contain¬ 
ing  different  chemical  potentials  of  oxygen,  such  as  SOFC 

m. 

7.5.  Measurement  of  reduced  negative  electrochemical 
potential  of  electrons,  cp,  by  embedded  reference  electrodes 

In  the  case  of  oxygen  ion  conductors,  the  probe  should 
exchange  only  electrons  and  not  oxygen  in  order  to  measure 
cp.  Since  the  ambient  atmosphere  contains  oxygen,  it  is  neces¬ 
sary  that  the  probe  be  completely  isolated  from  the  ambient. 
Fig.  3  shows  a  schematic  of  the  experimental  arrangement  for 
the  measurement  of  ip  in  a  solid-state  oxygen  ion  conducting 
electrolyte  by  embedding  a  platinum  reference  electrode.  The 
exposed  part  of  the  probe  is  sealed  off  from  the  atmosphere 
using  a  glass  which  is  impermeable  to  oxygen.  Voltage  is  mea¬ 
sured  between  a  probe  connected  to  one  of  the  electrodes  and 
the  embedded  probe.  When  the  probes  are  equilibrated  with  the 
sample, 


£e(Pt#l)  =  /Xg  °(electrodeI) 

(23) 

/xe(Pt#2)  =  pf. =x  (embedded  Pt) 

(24) 

The  voltage  measured  between  the  embedded  probe  and  elec¬ 
trode  I  is  given  by  [1]: 


AV=  Ae(Pt#2)  [  /Xe(Pt#l) 


fix~x  (embedded  Pt) 


+ 


/Xg_0  (electrode  I) 


(p(x)  —  cp( 0)  =  A  cp 


(25) 


Thus,  using  platinum  wires  and  platinum  reference  electrodes 
embedded  in  the  electrolyte,  the  reduced  negative  electrochem¬ 
ical  potential  of  electrons  (or  reduced  electrochemical  potential 
of  holes),  (p(x),  can  be  measured.  An  important  requirement  is 
that  the  embedded  probe  should  exchange  only  electrons  and 
not  oxygen. 


1.6.  Estimation  of  oxygen  chemical  potential  within  the 
electrolyte  from  the  measurement  ofcp 


Based  on  the  criterion  of  local  equilibrium,  we  have 
1  1 

Vcp  =  ——  V/i02-  +  —  V/*o2  (26) 

2e  \e 

At  open  circuit,  7i  +  7e  =  0.  If  the  electrolyte  is  a  predominantly 
oxygen  ion  conductor  (such  as  YSZ),  at  open  circuit,  V/xQ 2-  = 
0.  Eq.  (26)  then  reduces  to 

V<P  =  2  Vp IQ  (27) 

Ae 

Thus,  gradient  in  /x o2  can  be  readily  obtained  from  the  mea¬ 
surement  of  ip.  At  open  circuit,  from  Eq.  (27),  the  chemical 
potentials  of  oxygen  just  inside  the  electrolyte  at  the  cathode 
and  anode  are  given  by 

l4h  =  Mo2  -  4 e(<pl  ~  <fc)  =  Mo2  -  4 er\lt 

=  Mo2+4ere7i  (28> 

Mo2  =  /42  +  4^y‘  - /  )  =  fl g2  +  4<7e 

=  M02  -  4ere7i  (29) 


The  corresponding  oxygen  partial  pressures  within  the 
electrolyte  just  under  the  cathode  and  the  anode  are  given, 


Fig.  3.  A  schematic  showing  an  embedded  reference  electrode  in  a  solid  electrolyte  membrane.  The  reference  electrode  (probe)  is  sealed  off  from  the  atmosphere 
using  an  oxygen  impermeable  glass. 
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respectively  by 

C  ~  I  (  4e(^r  -  ^C) 

po2  =  po2^p[ - ^ — 


a  ~  n  {4e(<pa-  <pn) 
Po2  =  Po2  exP 


kBT 


Po2  exP 


Po2  exP 


\er%h 
'  kBT 


(30) 


(  4  erfIQ\ 

( tr) ,3i) 


The  corresponding  oxygen  partial  pressures  within  the  elec¬ 
trolyte  just  under  the  electrodes  are  given  by 


P(h  =  Po2  exp 


Po2  exP  l  - 


\e(r\  +  rf)(cpl  -  (pc) 
rfkBT 

kBT 


(39) 


In  YSZ-based  membranes,  the  electronic  leakage  is  neg¬ 
ligible.  Thus,  the  net  transport  of  oxygen  from  the  cathode 
to  the  anode  at  open  circuit  is  negligible,  and  Eqs.  (28)— (3 1) 
can  be  used  to  estimate  chemical  potentials  and  partial  pres¬ 
sures.  Abrupt  changes  in  both  are  expected  at  the  two 
electrolyte/electrode  interfaces,  the  magnitudes  of  these  changes 
being  dependent  upon  the  relative  values  of  transport  parameters 

m. 

If  the  electrolyte  exhibits  significant  electronic  conduction 
(such  as  rare  earth  oxide-doped  ceria),  then  V/xQ 2-  /  0  even  at 
open  circuit  and  significant  oxygen  permeation  may  occur.  In 
most  ceria  electrolyte-based  SOFC  operating  between  ~600  °C 
and  800  °C,  the  ionic  transference  number  of  the  cell  is  typi¬ 
cally  ~0.8.  The  transference  number  of  the  membrane,  including 
interfaces,  is  then  given  by 

t  —  £m  _  _ re  +  re‘  +  re _  _  _ 

E  rf  +  rf  +  rf  +  r  g  +  rf  +  rf  R\  +  Rq 

where 


Ri  =  rf  +  rf  +  rf  (33) 

Re  =  +  rf  +  ra  (34) 

and  Em  is  the  open  circuit  voltage. 

At  open  circuit,  the  internal  leakage  current  is  given  by 


E 

/leakage  “ 


=  -h  =  h 


(35) 


As  the  internal  leakage  current  is  not  negligible  in  ceria,  it  also 
means  that  V/xQ 2-  /  0.  The  knowledge  of  V/xQ 2-  is  required 
to  estimate  the  chemical  potentials  within  the  electrolyte.  Esti¬ 
mation  of  V/x02-  requires  the  knowledge  of  rf,  rf1,  and  rf .  The 
average  V /xQ 2-  (including  interfaces)  can  be  estimated  from  a 
cell  performance  test.  It  can  be  shown  that  the  voltage  vs.  current 
density  relation  is  given  by 


ER&  R[Re 

Ri  +  Rq  R[  +  Rq 


(36) 


Thus,  from  the  slope  and  the  intercept  of  this  curve,  both  R\ 
and  Rq  can  be  obtained.  The  rf 1  for  the  electrolyte  can  be  obtained 
using  separate  measurements  involving  electron  blocking  elec¬ 
trodes.  It  can  be  shown  that  the  oxygen  chemical  potentials  just 
inside  the  electrolyte  at  OCV  are  given  by  [7]: 

Mo2  =  m o2  -  4e(re  +  rf)4  =  p o2  +  4e(rt  +  riVi  (37) 

Mo2  =  <  +  4e(re  +  rf)/e  =  Mo2  _  4e(rl  +  ril)/i  (38) 


PO2  =  PO2  eXP 


=  Po2  exp 


/  4 e(rf  +  rg)(ya  - 

V  r*kBT 

(Mrf  +  r^lA 

V  kBT  ) 


(40) 


The  preceding  shows  that  the  knowledge  of  rf ,  rf ,  rf,  and  rf 
is  necessary  for  estimating  /Xq2  and  /Xq2  . 

2.  Experimental  procedure 

2.7.  Cell  fabrication 

Anode- supported  cells  were  fabricated  using  the  follow¬ 
ing  procedure.  Powder  mixture  containing  70wt.%  NiO  and 
30  wt.%  YSZ  (composition:  8  mol%  Y2O3-92  mol%  ZrC>2)  was 
made  to  which  some  carbon  black  was  added  to  create  poros¬ 
ity  upon  bum  off.  Additional  porosity  was  created  later  in 
the  anode  upon  reduction  of  NiO  to  Ni  during  cell  testing. 
Discs  of  3.3  cm  in  diameter  and  1.5  mm  in  thickness  were  die- 
pressed  and  fired  at  950  °C  in  air  for  1  h.  These  discs  were 
used  as  anode  supports.  A  slurry  containing  50  wt.%  NiO  and 
50  wt.%  GDC  (composition:  20mol%  Gd2O3-80mol%  Ce02) 
was  prepared  in  butyl  alcohol.  Anode  functional  layer  was 
applied  using  the  NiO  +  GDC  slurry  on  one  side  of  each  anode 
support  disc  by  drop  coating.  The  discs  were  again  fired  at 
950  °C  for  1  h.  Then,  the  first  electrolyte  layer  was  drop  coated 
using  a  slurry  of  GDC  followed  by  firing  at  1200  °C.  Platinum 
paste  was  applied  in  the  form  of  a  thin  and  a  very  narrow 
strip  over  the  electrolyte  layer  to  form  a  reference  electrode 
(Ref#2).  The  sample  was  heated  to  950  °C  to  remove  organics 
from  the  Pt  paste.  A  second  layer  of  GDC  was  drop  coated 
to  cover  most  of  the  surface  except  for  a  small  part  of  the 
Pt  strip,  which  was  left  exposed.  The  cell  was  again  fired  at 
1200  °C.  Another  reference  electrode  (Ref#l)  was  similarly 
applied,  followed  by  heating  to  950  °C.  Finally,  one  more  layer 
of  GDC  was  drop  coated  over  most  of  the  surface  leaving 
small  parts  of  both  electrodes  exposed.  The  cell  was  sintered 
at  1450  °C  for  10  h.  The  total  thickness  of  the  GDC  elec¬ 
trolyte  was  ~45  jxm  with  the  two  reference  electrodes  positioned 
~15  [Jim  from  the  cathode/electrolyte  and  the  anode/electrolyte 
interfaces,  respectively.  One  cell  was  made  with  a  total  of 
three  embedded  reference  electrodes;  one  (Ref#3)  ~15  pan 
from  anode/electrolyte  interface,  second  (Ref#2)  ~30  [xm  from 
anode/electrolyte  interface,  and  third  (Ref#l)  ~38  (xm  from  the 
anode/electrolyte  interface.  After  sintering,  cathode  functional 
layer  containing  a  mixture  of  LSC  and  GDC  was  applied  fol¬ 
lowed  by  firing  at  1100°C.  A  layer  of  LSC  was  then  applied 
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Glass 

Sealing 

REF  #1 


0 

0 

0 

REF  #2 
Glass  Sealing 


Fig.  4.  A  schematic  of  an  anode-supported  cell  with  embedded  reference  elec¬ 
trodes.  The  electrolyte  thickness  is  shown  exaggerated  with  respect  to  the  anode 
support.  AS,  anode  support;  AF,  anode  functional  layer  or  interlayer;  CF,  cathode 
functional  layer  or  interlayer;  CC,  cathode  current  collector. 


as  a  cathode  current  collector  followed  by  again  firing  at 
1100°C.  Pt  wires  were  attached  to  the  exposed  parts  of  the 
embedded  Pt  probes.  Exposed  parts  of  the  Pt  reference  elec¬ 
trodes  were  coated  with  a  paste  of  sealing  glass  in  ethylene 
glycol.  The  cell  was  fired  at  850  °C  to  fuse  the  glass  and  iso¬ 
late  the  exposed  parts  of  the  reference  electrodes  from  air. 
This  completes  the  cell  fabrication  procedure.  Fig.  4  shows 
a  schematic  of  the  cell  along  with  the  position  of  reference 
electrodes.  One  cell  was  made  with  two  embedded  probes  and 
a  thin  layer  of  YSZ  on  top  of  GDC  to  essentially  block  off 
electronic  current.  This  cell  was  used  to  measure  the  Nemst 
voltage. 

2.2.  Cell  testing 

Cells  with  embedded  Pt  probes  were  tested  using  a  specially 
designed  fixture.  The  cell  was  secured  at  the  end  of  an  alumina 
tube  with  the  anode  facing  the  inside  of  the  tube.  Mica  gasket 
was  used  as  a  glass-free  seal.  The  cell  was  held  in  position  by 
spring  loading  and  was  placed  inside  a  tube  furnace  with  the 
spring  being  outside  the  hot  zone.  The  cell  was  then  heated  to 
650  °C.  A  gas  mixture  containing  10%  H2  and  balance  nitrogen 
was  circulated  through  the  alumina  tube  over  the  anode  surface 
to  reduce  NiO  to  Ni.  Subsequently  the  anode  gas  was  changed 
to  —100%  H2.  Voltages  between  the  anode  and  the  reference 
electrodes,  and  between  the  anode  and  the  cathode  were  mea¬ 
sured  under  open  circuit  conditions.  Voltage  vs.  current  density 
polarization  curve  was  measured  on  the  cell  with  three  embed¬ 
ded  reference  electrodes.  OCV  was  measured  on  a  cell  with  a 
thin  YSZ  coating  over  GDC. 


2.3.  Measurement  of  the  conductivity  of  sealing  glass 

Three  discs  of  YSZ  were  made  by  die-pressing  followed  by 
sintering.  The  surfaces  of  the  discs  were  flat  ground.  A  sand¬ 
wich  was  formed  using  two  of  the  discs  with  a  layer  of  glass  in 
between.  The  sample  assembly  was  heated  to  850  °C  to  fuse  the 
glass.  Platinum  electrodes  were  applied  on  the  surfaces  of  the 
sandwich  sample  as  well  as  the  remaining  single  disc.  Two  probe 
dc  resistance  measurements  were  made  on  the  samples  at  650  °C. 
From  the  difference  of  the  two  measurements,  approximate  value 
of  the  ionic  conductivity  of  the  glass  was  obtained. 

2.4.  Measurement  of  total  conductivity  of  GDC  in  oxidizing 
and  reducing  atmospheres 

A  bar-shaped  sample  of  GDC  was  fabricated  by  die-pressing 
a  powder  of  GDC  and  sintering  at  1450  °C.  Platinum  strips  were 
painted  on  the  bar  in  a  four  probe  configuration.  The  sample  was 
heated  to  650  °C  and  its  four  probe  dc  resistance  was  measured. 
At  650° C  in  air,  the  predominant  transport  is  known  to  be  by 
oxygen  vacancies.  Subsequently  a  gas  mixture  containing  10% 
H2  +  balance  N2  was  introduced,  which  was  later  switched  to 
—100%  H2.  The  total  resistance  was  measured  as  a  function 
of  time.  Stable  values  in  air  and  hydrogen  were  recorded  as 
the  resistances  corresponding  to  the  two  atmospheres.  From  the 
sample  dimensions,  the  corresponding  resistivities  were  deter¬ 
mined. 

3.  Results  and  discussion 

3.1.  Micro  structures  of  cells  with  embedded  electrodes 

Fig.  5(a)  and  (b)  are  scanning  electron  micrographs  (SEM)  of 
a  tested  cell  showing  ESC  +  GDC  cathode  interlayer,  the  GDC 
electrolyte,  and  Ni  +  GDC  anode  interlayer.  Also  seen  in  the 
micrographs  are  the  two  embedded  Pt  electrodes,  Reference#  1 
(close  to  the  cathode)  and  reference#2  (close  to  the  anode).  The 
GDC  electrolyte  thickness  is  —45  pm.  Reference#  1  and  refer- 
ence#2  are  —  15  pm  from  the  cathode/electrolyte  interface  and 
—15  pm  from  the  anode/electrolyte  interface,  respectively.  Note 
that  the  GDC  electrolyte  appears  to  be  fully  dense  with  negligible 
porosity. 

3.2.  Measurement  of  the  conductivity  of  the  glass 

Two  probe  dc  resistance  measurements  were  made  at  650  °C 
on  a  YSZ  disc  with  Pt  electrodes  and  on  the  YSZ-glass-YSZ 
sandwich  sample  with  Pt  electrodes  (0.785  cm2  area).  The 
resistance  of  an  YSZ  disc  of  1.5-mm  thickness  was  mea¬ 
sured  to  be  —50  £2.  The  resistance  of  the  sandwich  sample 
was  measured  to  be  — 1 1,000  £2.  From  these  two  values, 
the  resistance  of  the  glass  was  estimated  as  —  1 0,900  £2. 
The  thickness  of  the  glass  layer  was  0.5  mm.  From  these 
results,  the  oxygen  ion  resistivity  of  the  glass  is  estimated 
to  be  -1.71  x  105Q  cm.  That  is,  the  ionic  resistivity  of  the 
glass  is  orders  of  magnitude  larger  than  that  of  YSZ  and 
even  larger  compared  to  GDC.  This  ensures  that  glass  seal- 
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Fig.  5.  (a)  A  scanning  electron  micrograph  (SEM)  of  the  cross-section  of  a 
cell  showing  the  anode  interlayer,  the  GDC  electrolyte,  the  cathode  interlayer, 
and  the  embedded  electrode  (reference  #1)  in  the  GDC  electrolyte  ~15  |xm 
from  the  cathode  interlay er/GDC  electrolyte  interface,  (b)  A  scanning  electron 
micrograph  of  the  cross-section  of  the  cell  showing  the  anode  interlayer,  the 
GDC  electrolyte,  the  cathode  interlayer,  and  the  embedded  electrode  (reference 
#2)  in  the  GDC  electrolyte  ~15  jxm  from  the  anode  interlayer/GDC  electrolyte 
interface. 


3.4.  Measurement  of  the  reduced  negative  electrochemical 
potential  of  electrons,  (p 

Since  GDC  is  a  mixed  ionic  electronic  conductor,  the  embed¬ 
ded  Pt  probe  can  exchange  both  oxygen  ions  and  electrons, 
or  effectively  exchange  oxygen.  However,  because  of  the  high 
ionic  resistivity  of  the  glass  used  for  sealing,  exchange  of  oxy¬ 
gen  is  suppressed.  Thus,  the  embedded  probes  should  ideally 
exchange  only  electrons  which  facilitates  the  measurement  of 
ip.  Fig.  6  shows  the  measured  voltages  between  the  anode  and 
Ref#2,  the  anode  and  Ref#l,  and  the  anode  and  the  cathode 
as  a  function  of  time  at  650  °C  with  100%  H2  as  the  anode 
gas  circulated  at  a  flow  rate  of  300  ml  min-1.  The  data  shown 
in  Fig.  6  are  after  the  establishment  of  steady  state.  The  OCV 
(anode/cathode)  is  ~0.85V.  The  anode  gas  was  not  bubbled 
through  water.  However,  some  H2O  was  formed  at  the  anode 
due  to  a  leakage  current  through  the  GDC  electrolyte  as  well 
as  possible  leaks  through  the  mica  edge  seal,  and  the  corre¬ 
sponding  po2  was  established  at  the  anode.  Fig.  6  shows  that 
the  voltages  were  very  stable  over  the  duration  of  the  test.  The 
voltages  between  the  anode  and  the  two  reference  electrodes 
were  very  small  (a  few  mV).  In  order  to  determine  if  this  was 
due  to  possible  shorting  of  the  probes  with  the  anode,  test  was 
also  conducted  with  10%  H2  +N2  as  the  anode  gas  circulated 
at  a  flow  rate  of  300  ml  min-1 .  Fig.  7  shows  the  results  of  these 
measurements.  In  this  experiment,  the  voltages  were  not  sta¬ 
ble  and  this  suggests  that  steady  state  had  not  been  established. 
The  OCV  was  between  ~0.7  V  and  ~0.8  V.  Note  that  voltages 
between  the  anode  and  Ref#2  and  the  anode  and  Ref#l  are 
nonzero — they  are  between  about  0.6  V  and  0.7  V.  This  obser¬ 
vation  suggests  that  the  reference  electrodes  were  not  shorted 
with  the  anode.  The  reason  for  some  drop  in  OCV  for  times 
>400  min  is  not  clear.  Also,  the  OCV  should  be  greater  than 
anode/reference  electrode  voltage,  but  part  of  the  time  it  was 
lower.  These  discrepancies  may  be  the  result  of  non- steady  - 
state  conditions  or  due  to  possible  incomplete  sealing  of  the 
probes  from  the  atmosphere.  After  ^700  min,  the  anode  gas 


ing  should  virtually  shut  down  any  oxygen  exchange  at  the 
probes. 


3.3.  Measurement  of  the  conductivity  of  GDC  in  air  and  in 
a  hydrogen  containing  atmosphere 

The  resistivity  of  the  GDC  bar-shaped  sample  at  650  °C  in 
air  was  measured  to  ~19.5£2cm  in  air.  GDC  is  a  predomi¬ 
nantly  oxygen  ion  conductor  in  air  at  650  °C.  Thus,  the  measured 
conductivity  is  attributed  to  predominantly  oxygen  ion  conduc¬ 
tion.  When  exposed  to  hydrogen,  the  resistance  of  the  sample 
rapidly  decreased  and  eventually  reached  a  stable  value.  The 
resistivity  of  the  sample  was  measured  to  be  ~4.55  Q  cm.  The 
oxygen  ion  resistivity  is  expected  to  be  essentially  independent 
of  oxygen  partial  pressure  and  thus  should  be  about  the  same 
as  ~  19.5  £2  cm.  From  these  two  measurements,  the  electronic 
resistivity  of  GDC  at  650  °C  in  hydrogen  is  estimated  to  be 
~5.935  £2  cm. 


0.8  -  Anode/Cathode 

0.7  - 


0.6  - 


0.2  - 


805; 

—  ■  —  Anode/Ref#2 

—  •  —  Anode/Ref#1 

—  *  —  Anode/Cathode 

O)  0.4  - 

TO 

£  0.3- 

0.1-  Anode/Ref#2  ~  Anode/Ref#1 

0-0  -  I  1  I  I  )  I  I  I 

300  325  350  375  400 


time(min) 


Fig.  6.  Measured  voltages  between  the  anode  and  Ref#l  and  the  anode  and 
Ref#2,  and  between  the  anode  and  the  cathode  at  650  °C  under  open  circuit 
condition:  fuel,  hydrogen;  oxidant,  air.  Note  the  very  small  voltages  between 
the  anode  and  the  reference  electrodes. 
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Fig.  7.  Similar  measurements  on  a  cell  with  initially  10%  H2  +  balance  N2  as  the 
fuel.  Relatively  high  voltages  between  the  anode  and  the  embedded  reference 
electrodes  show  that  the  reference  electrodes  were  not  physically  shorted  with  the 
anode.  When  the  anode  gas  was  changed  to  ~  100%  H2 ,  the  anode/reference  elec¬ 
trode  voltages  sharply  dropped  to  low  values,  consistent  with  the  development 
of  significant  electronic  conduction  in  the  GDC  electrolyte. 
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Fig.  9.  Measurements  of  voltages  on  a  cell  with  three  embedded  reference  elec¬ 
trodes.  Measured  voltages  between  the  anode  and  Ref#l,  Ref#2,  and  Ref#3, 
and  between  the  anode  and  the  cathode  at  650  °C  under  open  circuit  condition: 
fuel,  hydrogen;  oxidant,  air.  Inset  shows  an  enlarged  vertical  scale.  Note  that  the 
voltage  between  the  anode  and  Ref#l  is  the  highest  (~30mV).  Slightly  higher 
voltage  between  the  anode  and  Ref#3  than  between  the  anode  and  Ref#2  is  likely 
experimental  scatter. 


was  switched  to  100%  H2.  Almost  instantaneously  the  volt¬ 
ages  between  the  anode  and  Ref#2  and  between  the  anode  and 
Ref#l  dropped  to  very  low  values  (a  few  millivolts).  The  voltage 
between  the  cathode  and  the  anode  (OCV),  by  contrast,  sharply 
increased  followed  by  a  slow  rise  to  about  0.85  V  in  about  a  cou¬ 
ple  of  hours.  The  observation  that  the  voltage  difference  (A cp) 
between  the  anode  and  the  reference  electrodes  is  very  small  is 
consistent  with  the  development  of  significant  electronic  con¬ 
duction  in  GDC.  Fig.  8  shows  the  variation  of  cp  through  the 
GDC  electrolyte. 

Fig.  9  shows  Acp  measurements  on  a  cell  with  three  embedded 
probes;  Ref#3  is  ~15  jx m  from  the  anode/electrolyte  inter¬ 
face,  Ref#2  is  ~30  |xm  from  the  anode/electrolyte  interface, 
and  Ref#l  is  ~38  |xm  from  the  anode/electrolyte  interface  or 
p,m  from  the  cathode/electrolyte  interface.  The  anode  gas 
was  ~100%  H2.  The  OCV  was  ~0.84  V.  The  measured  volt¬ 
ages  between  the  anode  and  Ref#2  and  between  the  anode  and 
Ref#3  were  very  small.  The  measured  voltage  between  the  anode 
and  Ref#l  was  also  small  (~30mV)  but  still  somewhat  larger 


than  between  the  anode  and  Ref#2  or  the  anode  and  Ref#3.  The 
measurements  show  that  cpu,  <pa,  <pref3,  and  cpVQf2  are  about  the 
same,  while  <prefl  is  slightly  larger.  In  order  to  determine  cpc,  it 
was  assumed  that  (p  varies  linearly  between  Ref#2  and  the  cath¬ 
ode/electrolyte  interface.  Since  Ref#l  is  symmetrically  located 
between  the  cathode/electrolyte  interface  and  Ref#2,  cpc  —  cpmfl 
should  be  approximately  about  the  same  as  cpr&fl  —  <pref2.  The 
measured  open  circuit  voltage  was  ~0.84  V  and  the  difference 
cpl  —  cpc  is  estimated  to  be  ~0.79  V.  By  contrast,  the  difference 
cpmi 3  —  cp11  was  measured  to  be  only  ~10  mV.  Thus,  for  all  prac¬ 
tical  purposes,  the  Acp  across  the  anode/electrolyte  interface 
(<pa  —  cp11)  is  negligible.  If  we  assume  cpa  —  cp11  ~  <pref3  —  cp11,  then 
this  drop  in  voltage  is  ~10  mV,  or  cpa  —  cp11  =  ra/e  =  10  mV.  By 
contrast,  at  the  cathode,  cp1  —  cpc  =  =  0.79  V.  These  results 

show  that  the  electronic  resistance  of  the  cathode/electrolyte 
interface  is  much  larger  than  that  of  the  anode/electrolyte 
interface.  Fig.  10  shows  a  plot  of  the  cp  as  a  function  of  posi¬ 
tion  in  the  GDC  electrolyte  of  the  cell  with  three  embedded 
probes. 
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Fig.  8.  A  plot  of  the  measured  cp  corresponding  to  the  data  given  in  Fig.  7  with 
~100%  H2  as  the  anode  gas. 
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H.-T.  Lim,  A.V.  Virkar  /  Journal  of  Power  Sources  180  (2008)  92-102 


101 


3.5.  Estimation  of  oxygen  chemical  potential  or  oxygen 
pressure  using  the  measurement  of  ip  at  open  circuit 

Voltage  vs.  current  density  test  was  conducted  on  the  cell  with 
three  embedded  electrodes.  The  measured  voltage  vs.  current 
density  plot  is  shown  in  Fig.  1 1 .  From  the  cell  test,  the  magnitude 
of  the  slope  is  identified  with  RiRJ(Ri+Rf)  and  is  measured  as 
~0.31  Q  cm2.  The  open  circuit  voltage  is  given  by  ERe/(R{  +Re) 
and  is  equal  to  0.84  V.  Cell  voltage  was  also  measured  on  a  GDC- 
based  cell  with  an  YSZ  coating  which  essentially  blocks  off  the 
electronic  current.  This  cell  voltage  was  measured  as  1.05  V 
at  650  °C,  which  is  assumed  to  be  the  Nernst  voltage,  E.  From 
these  measurements,  the  R{  and  Re  are  estimated  as  ~0.39  Q  cm2 
and  ~  1 .56  ^2  cm2,  respectively.  From  these  values,  the  leakage 
current  is  estimated  as  |/i|  =  |/e|  =0.54  A  cm-2.  Using  the  mea¬ 
sured  ionic  conductivity  of  GDC  at  650  °C  of  ~0.05  Scm-1, 
(Aq2—  —  Aq2-  )/2e  is  estimated  as  ~0.0474  V,  in  which  the  ionic 
conductivity  of  GDC  is  assumed  to  be  independent  of  po2 . 

In  what  follows,  we  will  assume  that  the  area  specific  ionic 
interface  resistance  at  the  anode/electrolyte  interface  is  negli¬ 
gible,  that  is  rf  ^  0.  The  observation  that  <pa  ~  ip11  also  means 
rf  ~  0.  In  such  a  case,  the  /Xq2  —  /Xq2  =  A/x q2  ~  0.  That  is, 
the  chemical  potential  of  oxygen  nearly  equilibrates  across  the 
anode/electrolyte  interface.  Indeed,  Raman  spectroscopic  stud¬ 
ies  by  Mineshige  et  al.  in  a  GDC  electrolyte  of  4-mm  thickness 
showed  that  the  chemical  potential  of  oxygen  in  the  electrolyte 
at  the  anode/electrolyte  interface  was  essentially  the  same  as 
in  the  fuel  [22].  This  observation  shows  that  indeed  rf  ~  0  is 
a  reasonable  assumption.  With  this  assumption,  and  the  mea¬ 
sured  ip  in  the  GDC  electrolyte,  the  variation  of  oxygen  chemical 
potential  within  the  GDC  electrolyte  can  be  readily  estimated. 
Fig.  12  shows  the  plots  of  ip  (measured  by  embedded  probes), 
/x02—  (estimated  from  the  ionic  conductivity  of  GDC,  the  leak¬ 
age  current  and  the  assumption  rf  ~  0),  and  In  po2  which  is 
a  measure  of  /xo2  (calculated  from  /x0 2-,  ip,  and  the  present 
results  showing  that  ra  ~  0).  The  estimated  value  of  the  oxygen 
partial  pressure  in  the  GDC  electrolyte  just  inside  the  cath¬ 
ode,  Pq2  ,  is  estimated  to  be  ~1.99  x  10-22atm,  while  that  at 
the  cathode  is  ~0.21  atm;  about  21  orders  of  magnitude  drop 


Fig.  1 1 .  A  performance  test  at  650  °C  on  the  cell  with  three  embedded  reference 
electrodes. 


REF  m 
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Fig.  12.  Plots  of  cp,  jl02-/2e,  and  In  po2  (or  effectively  /zo2)  through  the  GDC 
electrolyte.  Note  the  sharp  drop  in  po2  across  the  cathode/electrolyte  interface; 
~0.21  atm  at  the  cathode  and  ~1.99  x  10-22  atm  in  the  GDC  electrolyte  just 
across  the  cathode/electrolyte  interface. 

in  oxygen  partial  pressure  just  across  the  cathode/electrolyte 
interface.  This  is  quite  a  remarkable  result  in  that  it  shows  that 
there  is  an  abrupt  change  in  /x o2  across  the  cathode/electrolyte 
interface;  that  is  A /Xq2  =  /Xq2  —  /Xq2  /  0,  and  in  fact  is  rather 
large.  A  linear  extrapolation  of  In  po2  between  ref  2  and  ref 
1  up  to  the  cathode/electrolyte  interface,  which  is  an  approx¬ 
imation,  was  assumed.  Even  if  some  nonlinearity  is  allowed 
for,  it  is  still  clear  that  there  is  an  abrupt  change  in  po2  across 
the  cathode/electrolyte  interface.  These  results  also  show  that 
the  often  made  assumption  of  the  equilibration  of  chemical 
potential  of  oxygen  across  a  gas/solid  interface  when  there  is 
a  net,  nonzero  A/xo2  across  the  membrane,  is  generally  not 
true.  We  have  also  seen  that  A cpc  =  (pl  —  <pc  ~  0.79  V.  This  means 
the  electronic  area-specific  resistance  at  the  cathode/electrolyte 
interface  is  given  by  rf  =  A<pc/|/e|  or  rf  ~  0.79/0.54  £2  cm2  or 
~  1.463  £2  cm2. 

Now  the  Ri~0.39£2cm2  and  ra~0£2cm2.  From  the 
measurement  of  conductivity  of  GDC  in  air,  the  estimated 
value  of  the  electrolyte  ionic  area- specific  resistance  is  rf  ~ 
0.09  Q  cm2.  From  these  measurements,  since  R{  =  rf  +  rf  + 
rf,  the  area- specific  cathode/electrolyte  ionic  interface  resis¬ 
tance  is  estimated  as  rf  ~  0.3  Q  cm2.  The  performance  of  GDC 
as  a  good  electrolyte  for  SOFC  applications  is  thus  attributed 
to  the  large  value  of  rf,  electronic  resistance  at  the  cath¬ 
ode/electrolyte  interface,  which  leads  to  moderately  high  ionic 
transference  number  for  the  cell.  The  present  results  also  thus 
suggest  that  virtually  all  of  the  GDC  electrolyte  thickness  even  in 
a  thin  electrolyte  anode- supported  SOFC  is  a  mixed  ionic  elec¬ 
tronic  conductor  and  much  of  the  voltage  drop,  A <p,  and  much 
of  the  chemical  potential  of  oxygen  drop,  A/xo2 ,  drop  across  the 
cathode/electrolyte  interface. 

In  the  work  by  Mineshige  et  al.,  analysis  of  their  data  was 
made  under  the  assumption  that  the  /xo2  equilibrated  across 
both  cathode/electrolyte  and  anode/electrolyte  interfaces;  that 
is  /Xq2  ~  /Xq2  and  /Xq2  ~  /Xq2  [22] .  Their  experimental  results 
on  the  anode  side  were  in  agreement  with  the  assumption 
/Xq2  ~  /xg2.  However,  examination  of  their  Fig.  5  shows  that 
/Xq2  /  /Xq2  and  that  /Xq2  <  /Xq2,  and  thus  not  in  agreement 
with  the  assumption  /Xq2  ~  /Xq2  made  in  their  work.  That  is, 
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in  their  work  also  there  was  a  sharp  drop  in  /xo2  across  the  cath¬ 
ode/electrolyte  interface,  which  is  not  in  agreement  with  their 
(and  commonly  made)  assumption  of  the  equilibration  of  chem¬ 
ical  potential  across  the  gas/solid  interface.  However,  the  abrupt 
change  in  fio2  observed  in  their  studies  and  in  the  present  work  is 
in  complete  accord  with  the  present  work  and  analysis  presented 
previously  [7]. 

4.  Summary 

Anode- supported  cells  with  GDC  electrolyte  of  thickness 
~45  p,m,  LSC  +  GDC  cathode  interlayer,  Ni  +  GDC  anode  inter¬ 
layer,  and  Ni  +  YSZ  anode  support  were  fabricated  with  Pt 
probes  embedded  within  the  GDC  electrolyte.  Up  to  three  Pt 
probes  were  embedded  within  the  electrolyte.  The  probes  were 
sealed  off  from  the  atmosphere  using  a  glass  impermeable  to 
oxygen.  Voltages,  A cp,  between  the  anode  and  the  embedded 
probes  and  between  the  anode  and  the  cathode  were  measured  at 
650  °C  with  air  circulated  past  the  cathode  and  hydrogen  past  the 
anode.  It  was  observed  that  there  was  negligible  voltage  differ¬ 
ence  between  the  anode  and  the  embedded  reference  electrodes, 
and  much  of  the  voltage  dropped  across  the  cathode/electrolyte 
interface.  This  drop  is  attributed  to  a  large  electronic  area- 
specific  resistance,  r£,  across  the  cathode/electrolyte  interface. 
Based  on  these  results,  the  chemical  potential  of  oxygen,  /io2 , 
through  the  GDC  electrolyte  was  mapped  out  as  a  function  of 
position.  The  results  show  that  much  of  the  chemical  poten¬ 
tial  of  oxygen,  A/xo2,  also  drops  across  the  cathode/electrolyte 
interface.  The  present  results  show  that  in  an  SOFC  made  with 
a  mixed  ionic  electronic  conducting  ceria  electrolyte,  such  as 
GDC,  only  an  extremely  thin  layer  at  the  cathode/electrolyte 
interface  functions  as  an  electrolyte. 
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